Comparative study of surface plasmon scattering by shallow ridges and grooves 
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We revisit the scattering of surface plasmons by shallow surface defects for both protrusions and 
indentations of various lengths, which are deemed infinite in one-dimension parallel to the surface. 
Subwavelength protrusions and indentations of equal shape present different scattering coefficients 
when their height and width are comparable. In this case, a protrusion scatters plasmons like 
a vertical point-dipole on a plane, while an indentation scatters like a horizontal point-dipole on 
a plane. We corroborate that long and shallow asymmetrically-shaped surface defects have very 
similar scattering, as already found with approximate methods. In the transition from short shallow 
scatterers to long shallow scatterers the radiation can be understood in terms of interference between 
a vertical and a horizontal dipole. The results attained numerically are exact and accounted for 
with analytical models. 
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I. INTRODUCTION 

Surface Plasmon Polaritons (SPP) are electromagnetic 
bound modes responsible for the transport of light at 
the interface separating a metal from a dielectric. Their 
ability to confine light at an air-dielectric interface of- 
fers the prospect of developing a new technology con- 
sisting of photonic nano- devices 1-4 . Active research is 
currently focusing on the possibility of achieving control 
over the propagation of SPPs by means of optical ele- 
ments that would couple or decouple light to them 5-9 . In 
order to conceive optical elements (lenses, mirrors, beam- 
splitters) able to manipulate SPPs propagation, we need 
to learn more about the interaction of surface plasmons 
with a sub-wavelength modification of the underling di- 
electric metal interface. Indeed the interaction of SPPs 
with surface sub- wavelength defects on a metal surface is 
of great interest from a theoretical standpoint. 10 ' 11 
In this article we shall study scattering of SPPs by a 
shallow surface defect. We will consider both indenta- 
tions of the metal surface (grooves) and protrusions on it 
(ridges). We shall only deal with bi-dimensional defects, 
which are deemed infinite in one dimension parallel to 
the interface (the y-direction). Different aspects of this 
problem have been studied before with a variety of nu- 
merical techniques 12-21 . Here we present a systematic 
comparison between the different scattering coefficients 
and provide both analytical expressions and qualitative 
explanations. 

It must be noted that in a previous work we presented 
such a comparison 22 but within an approximate numer- 
ical scheme. Within that framework it was found that 
ridges and grooves exhibited the same scattering, when- 
ever they are shallow enough. Here we will revise that 
result, which turns out to be valid only for long (elon- 
gated) defects. The mistaken outcome of Ref. [22] for 
short defects may be traced back to the breakdown of 
the assumption of small curvature in the defect geometry 
that was made there. In this paper we solve the Maxwell 
equations through a discretization method, which does 



not assume the previous approximation and whose accu- 
racy depends only on the discretization mesh. We found 
that, as in the previous work, long asymmetric ridges or 
grooves with the width much larger than the depth, do 
scatter very similarly. However square shallow defects 
manifest a different scattering efficiency, differing in the 
relative radiative loss and radiation pattern. The lack of 
distinction between these two cases did not emerge in the 
previous approximate treatment. On the whole the prob- 
lem needs to be revisited so as to: i) substantiate why 
the approximate result does work in the case of elongated 
defects, ii) point out what is the correct result in the case 
of shallow and short symmetric defects, and in) explain 
qualitatively how the scattering properties of short and 
shallow symmetric defects are gradually transformed into 
the scattering properties of elongated defects, as the as- 
pect ratio of the defect increases. 

This paper is organized as follows. In Sec. II we state the 
basic assumptions on the scattering system as well as the 
solution method. In Sec. Ill we rearrange the asymptotic 
expansions of the far-field to produce the scattering coef- 
ficients. Namely we express the far-field and the related 
Poynting vector in terms of the field inside the defect. 
Still in this section we look at an approximation for the 
scattering coefficients of shallow ridges. In Sec. IV we 
explain that, in general, we cannot quantitatively repre- 
sent a scatterer (however small) by one mesh. We explain 
how we associate a small symmetric ridge or groove to 
a point dipole. In Sec. V we look at exact numerical 
results for the scattering of shallow defects of various 
horizontal lengths. We analyze these results and, in the 
case of square defects, we associate a ridge to a vertical 
dipole and a groove to a horizontal dipole. In Sec. VI 
we produce an analytical model that explains the radia- 
tion pattern of the surface plasmons scattered by small 
square ridges and grooves. In Sec. VII we look at the 
solutions for the case of shallow and long defects and 
we present a clear-cut interpretation to support the re- 
sults of the previous treatment 22 . Finally in Sec. VIII 
we explain qualitatively that the aspect ratio of the de- 
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feet determines the orientation of the held induced in a 
shallow defect. 



II. THE SCATTERING SYSTEMS 
CONSIDERED 

The considered defects are infinite in the y-dimension 
and shallow with depth h << A, where A is the free 
space wavelength. The defects are going to be illumi- 
nated by a monochromatic surface plasmon at normal 
incidence e spp , associated to an impinging energy flux 
S spp , defined and derived in Appendix A. Therefore, 
only radiation into p-polarized(TM) waves needs to be 
considered. After we drop, out of symmetry, the y- 
dependence on the whole problem the field is expressed 
as: E(r,t) = E(x, z)e~ tult . The wavevector in vacuum is: 
g = 27r/ A , where uj = eg. The material making the slab 
shall be lossless silver 23 , that is: e = 3?{e J 4 9 (A)}. Absorp- 
tion is neglected as we consider non-resonant defects with 
widths much smaller than the SPP propagation length. 
As represented in Fig. 1, we shall be expressing the source 
orientation in a cartesian basis (u x ,u z ), and the scat- 
tered fields in a right-handed orthogonal polar basis: 

Ufj — cos au x + sin au z (1) 
u Q = — sin au x + cos au z (2) 

Finally a question of notation: throughout we shall re- 
fer to a bi-dimensional point-source simply as a dipole, 
but it is meant that their emission as all of the fields 
are cyclical in the y-direction. As represented in Fig. 1, 
each object lying in the vacuum semi-space shall be la- 
beled by the superscript v = 1 while any object lying 
in the metal shall be labeled by the superscript v = 2. 
In particular, scattering quantities related to ridges have 
the superscript v = 1 while the ones related to grooves 
have the superscript v = 2. The field within the cross- 
sectional area of the ridge is labeled E r (r') and the one 
within that of the groove is labeled E s (r') . 



III. SCATTERING COEFFICIENTS 

The Green tensor approach is a standard method to 
solve electromagnetic scattering problems 1 1,24 ~ 30 . Our 
first task in this section, is to arrive at an explicit ex- 
pression for the scattered electric far-field. This is at- 
tained by propagating the field induced by a dipole den- 
sity PW(r') = Ae E r (r') (where Ae = e - 1) inside 
the area of a ridge, to a point R very far from the 
source. For a groove we have the same relation be- 
tween polarization and field (except for a change of sign) 
p( 2 )(r') = -Ae E 9 (r'). To propagate the field from any 
of the two, we use the standard formula 11 : 

E S (R)= 5 2 / dr'G(R,r')-P(r'). (3) 

J A 




FIG. 1: Schematic representation of the scattering systems 
considered. A ridge is a one dimensional defect located in air 
and labeled by the index v = 1. A groove is a one dimensional 
defect located in the metal and labeled by the index v — 2 



Where G(R, r') is the Green tensor for the air- metal 
background. The Green tensor propagates the emission 
of a point source at r' to the distant point of the R. One 
of the advantages of the Green tensor technique is that 
once the fields inside the defects E(r') (and thus P(r')) 
are computed numerically the asymptotic expansions of 
scattered fields become analytic. This takes us to our sec- 
ond task, which is making a direct connection between 
the orientation of the induced polarization inside the de- 
fects and the far-field radiation pattern, and in so doing 
define the scattering coefficients. 

First of all, finding the scattered electric far- field E S (R) 
requires the asymptotic expansions of the Green tensor. 
The derivation is sketched in the Appendix C. In what 
follows we give some simplifying rearrangements that will 
let us focus directly on the angular radiation pattern of 
surface defects. 



A. Scattering into Radiative Modes 

The asymptotic Green's tensor in the radiative zone for 
either a ridge or a groove can be written in a compact 
form as: 

e i{gR+ir/4) 

G^(R^ oo,a,r') = — e~ %!,x cosa x 

x e~ ih >' GM(a,r'). (4) 
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In such form we can factor the asymptotic scalar green 
function out of the dyadic part of the Green tensor. From 
eq.(3), the direction of E S (R) results from superposition 

of G { Z\a,r')-P v (r'), the emission from all induced point 
polarization elements, or dipole density elements. Yet the 
direction of each contribution (a, r') • P" (r') must be 
independent of r'. In other words since electromagnetic 
waves are transverse waves in vacuum, far from their 
source, the field emitted by a dipole G^)(a, r') • P iy (r') 
must be proportional to u a . In fact using the stan- 
dard asymptotic expansions (see the Appendix C) we can 
write: 

GW(«,r').F(r') = -[$M(«,r').F(r')] u a . (5) 
Where for a ridge: 

=k+(a)+k;( a )r p (a)e 2 ^' sinQ . (6) 
and for a groove: 

& 2 Ha) =t^ 2 Ha)k pm (a). (7) 

The vectors k p (a) are p- waves defined in vacuum, while 
k™ ± (a) are defined in the metal. A reminder of their 
expressions at normal incidence, in terms of the angle 
a of Fig. 1, is reported in the Appendix B, along with 
the expression for the Fresnel reflection and transmission 
coefficients: r p (a), t p 1,2 \a). 

We are now in a position to write the expressions for the 
radiative fields. Plugging eq.(4) and eq.(5) into eq.(3) 
we can separate the electric far field dependence into its 
radial and angular parts as: 



EM(fl,a) = - 



J(gR-*/4) 



E^(a)u a . (8) 



Here the angular amplitude can be written as: 

E^(a)=g 2 Q^ d (a) (9) 

where ©^(a) is the scattering coefficient into radiative- 
modes: 

e^(a) = J d r ' e -v'co SV .e ) ^M( Q;r '). P M( r ') 



(10) 



In the last expression the scattered field in the far zone 
consists of a cylindrical wave, transverse to the direction 
of propagation u#, and with a net angular amplitude de- 
termined by the integral over the source region Q^J d (a). 
The latter is actually the important bit in the formula 
as its squared module determines the radiation pattern. 
As seen from eq.(10) this angular amplitude results from 
the superposition of each scattering element taken with 
its own amplitude, phase and optical path in analogy to 
how an antenna array determines its effective radiation 



pattern. The radiation is given by the intensity or Poynt- 
ing vector in the far field. Accordingly the differential 
angular scattering cross-section is: 

d^M) \E^\R,a)?R _ g 3 . |2 

a - o - a \ KJ rad\ a )\ ' 

(11) 

Finally, the net radiative loss <J ra d is defined as the inte- 
grated angular radiation: 



&rad 



L 



180 ° da (v) Aa) 
oa 



(12) 



B. Shallow defects and Green's tensor boundary 
conditions 



Whenever the height of the defect is small enough, typ- 
ically much smaller than the wavelength of the incident 
light, we can make the approximation g\r'\ << 1. That 
allows for some simplification for the angular amplitude 
of a scattering element above the surface. Consider: 

*W(a,r') = (k+( a )+k;( a )r p ( a )e 2 ^' sinQ )^ 

~ (k+(a) + k^(a))r p (a) (13) 

Hence, for shallow defects the Green Tensor dependence 
of cq.(4) is entirely given by the exponential factors 
e -igx cosa e -ik z z ^ f QT j-, th a source in the vacuum semi- 
space and a source in the metal semi-space. Indeed this 
turns out to be a major simplification for the relative am- 
plitude of the scattering elements in the air semi-space, 
which we shall perform in detail SectionVI. 
Before that we need to highlight the relation between 
the Green tensor of a defect on the metal slab and in 
the metal slab, under this approximation. Such relation 
emerges from the boundary conditions for the Green's 
tensor at the interface, which are: 



G(R, x', z' = 0+) - G(R, x', z' = 0") 



G(R, x', z' = 0+) - e G(R, x', z = 0~) 



■ u x = 



u 2 



(14) 



(15) 



Notice that, in the unperturbed system, space is transla- 
tionally invariant in the horizontal direction x and this is 
reflected in is the x-component of the vector in eq.(13). 
Because of eq.(4) and eq.(5), we can turn eq.(14) into: 



(16) 



The presence of surface charges at the interface im- 
plies, from cq.(15), that the z-components of the vector 
&("\a) on either sides of the interface have the relation: 



(17) 
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C. Scattering into Surface Plasmons 

Let us derive the scattering coefficient into surface 
plasmon modes. Note that, in this one dimensional prob- 
lem, scattering will be into both the forward surface 
plasmon e spp+ (r), propagating in the positive x direc- 
tion and the backwards plasmon e spp _(r) propagating 
in the negative x direction, as defined in the Appendix 
A. The emission by a point dipole or a point polar- 
ization element must result into a plasmon final state: 
G P ±(R, r') • P(r') oc e spp ±, as shown in the derivation 
sketched in the Appendix C. The asymptotic Greens 
tensor for a source upon (1/ = 1) or in (1/ = 2) the metal 
is: 



G^(R,r')-P (l V) =^^x 
(ei(r'))*-PM(r')] e spp± (R). (18) 



Notice that ^e| pp± (r')^ complies with eq.(14) and 

eq.(15). Consequently the field of the scattered plasmons 
are: 

^p ~ 9 a ^p± ^ S PP± 

^ ^spp 

Opt = / dr'e* spp± (v').P(v') 

J 



(19) 



Furthermore the magnetic field related to the field scat- 
tered into SPPs is : 



H p ± = -- 



2S. 



e 



p± "spp± 



(20) 



spp 



where h spp is the magnetic field of a SPP, as proved in 
the Appendix A. 

Now, if the the source P(r') is produced by an incident 
surface plasmon field (as is our case), we can define the 
scattering cross-section of into SPPs as: 



E p± x H* ± • u x 
e p± x h* ± ■ u x 



fJ 



2S 



-e p± 



spp 



(21) 



Finally, we can define the total scattering cross-section, 
which in the lossless case is equivalent to the extinction 
cross-section: 



°~p + a p + °~rad 



(22) 



IV. 



RAYLEIGH-LIMIT: CAUTIONARY 
REMARKS 



Next we are going to develop solutions to point sources 
in a metal plane background. However one question may 
be raised : how do we associate the field induced by a sur- 
face plasmon inside a ridge or a groove to a point dipole? 



The answer is the argument of this section. 
When the field inside a defect is obtained by mesh dis- 
cretization we assume that the field inside a single mesh 
is uniform, and deviations from the field at its center are 
deemed negligible. Yet, in general, the field in a defect, 
cannot be represented by the field at its center alone. Let 
us explain a little bit further this point. 
For simplicity let us consider a defect in a homogenous 
medium with dielectric constant £{,, but the argument is 
the same in other backgrounds. As usual 31 , the field at 
every mesh is found by solving self-consistcntly a system 
of N coupled equations: 

E(r0 = E 6 (r0 + g 2 ]T G fc (r, - Tj ) ■ Ae Efo)-^ + 



+g 2 M • Ae E(i-j) AeE(r,) (23) 

£b 

where i = l,N and j = 1,N and E(r^) is the field at the 
mesh center. L is a term related to the depolarization of 
light and comes about from the quasi-static contribution 
of the Green tensor. M is a correction term to the Green 
tensor in the region of the scatterer useful to improve the 
accuracy of the calculation, when the inhomogeneity is 
discretized 11,32 . 

In practice, the number of mesh points N is increased 
until the calculation converges to the required precision. 
Then scale variations ~ y/A/N of E(r) are properly rep- 
resented in the solution. In the Rayleigh limit, for a 
defect of area A so small that g 2 A « 1, the scatterer 
behaves like a point source or a point dipole and the 
background field (in this case the illumination) can be 
considered uniform over A: Ef,(r) = E& . Exceptionally, 
for a circular defect in a homogenous medium with dielec- 
tric constant Sb, the net field at any point Vi converges 
to: 



E = E 5 AeE 

£b 



(24) 



This is because for the field inside an infinitesimal (very 
sub-wavelength) circular shape is actually uniform and 
thus scattering by such circular defects can be described 
by one mesh. In fact the extinction coefficient 33,34 can 
be derived from the field at the center alone: 



0~xtn — g^$ 



f dr'AeE*(r') • E(r') 

J A 



U A 

= A3 [AeE£ • E] 



(25) 
(26) 



To prove this numerically we have calculated a xtn for a 
cylinder represented by a single mesh, as in eq.(24), and 
illuminated by a plane wave. First of all we have checked 
that the one- mesh cross-section of eq.(26), coincides with 
the Mie theory result. Secondly, we have subdiscretized 
the cylinder into square meshes as rendered in the inset 
of Fig. 2. As also rendered in the figure, applying eq.(23) 
we found that, as the number of meshes grows, the scat- 
tering cross-section calculated by the collection of meshes 
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in this paper. We find that, for a defect above the sur- 
face in the optical range, the relative error is about 40%, 
while it can reach 50% for a defects below the surface. 
With very small non-elongated ridges and grooves, such 
that w/X ~ h/X << 1, the equivalent point dipoles are 
attained by averaging the fields over the area of the de- 
fects as follows: 



p (1) = AeE r A = Ae / dr'E r (r') 

J A 

= -AeE 9 A = -Ae / dr'E ff (r'K 

J A 



(28) 



(29) 



Accordingly if we set P(")(r) = 8(r — r')p'^ eq.(10) and 
eq.(19) for small non-elongated defects become: 



eg (a) = feg(0) 



P 



(30) 
(31) 



FIG. 2: Dependence of the extinction coefficient on the num- 
ber of meshes used in the calculation, for a square defect with 
w — h = lnm in vacuum, and a circular defect of the same 
area, illuminated by a plane wave. The dashed line represents 
extinction coefficient calculated with the Mie theory for the 
circle. The dielectric constant in the defect is e — —19.89 at 
the wavelength of 700nm. The inset represents the geometry 
of a discretized circle when inscribed in a square represented 
by 30x30 mesh points. 



eq.(25) converges to the initial value of one single mesh 
of eq.(26). However the field inside of a square scatterer 
can never be uniform if it is to satisfy real boundary con- 
ditions even in a homogenous medium or vacuum. Thus, 
it can not be faithfully described by one mesh. This is 
illustrated in Fig. 2, which renders the extinction coeffi- 
cient for a square defect of the same area as the circle. 
As it turns out, the converged value is ~ 27% larger than 
that obtained by the one-mesh approximation. Remark- 
ably this error is not reduced with the defect size: we 
obtained the same error for squares with side hnm or 
0.5nm. This is just for reference in the optical range, 
since we found that the error actually depends on type 
of defect and on the dielectric constant. 
However, even if the field is not uniform, a small defect in 
the Rayleigh limit can be represented by a point source 
at the center of the mesh, with its field equal to the av- 
erage field over the mesh E = (l/A) f. dr' E(r'). 
Indeed if the variation of Ef,(r) is negligible over the area 
of the defect we have: 

a xtn = AS [AeEl • E] (27) 

So the object behaves as a point-dipole p = AAe~E. 
The previous results were for a homogeneous background, 
but they also hold for the inhomogeneous one considered 



V. NUMERICAL RESULTS 

As an illustration consider a square ridge and a groove 
of side w = h = lOnm. We have calculated the scatter- 
ing into radiative modes and SPPs without associating 
the defect to a point dipole but rather using eq.(ll) and 
eq.(21). In this case the major task is computing the 
Green's tensor for the plane metal surface required to at- 
tain the exact field within the surface defect. This can 
be achieved following the prescriptions of Ref. [35,36]. 
Similar numerical results for the case of shallow grooves 
were found in Ref. [19] using a different computational 
technique. 

The out of plane radiation pattern of a surface plasmon 
scattered by such defects is given in Fig. 3. 
Calculations show that, for symmetric defects, the net 
radiative loss is greater for a groove than for a ridge. 
This is so because, while both the scattering into SPPs 
and the radiation close to the surface (at a = 0, 180°) 
are similar, their radiation patterns greatly differ normal 
to the surface (a = 90°), where the groove radiation is 
maximum while the ridge radiation goes to zero. The 
ridge radiation pattern is distributed into two lobes on 
either sides of a — 90° but the groove radiation pattern 
forms a single lobe. This is one of our main result and 
shall be analyzed in detail in the next section. The result 
is not in agreement with those obtained in the approx- 
imate treatment Ref. [22]. We associate the discrepancy 
to the breakdown of the condition that the curvature of 
a short and shallow defect does not vary rapidly, used in 
that work. 

Notice the fraction of energy scattered into SPPs, i.e 
of eq.(21), is large. The values of are represented 
by the horizontal lines of Fig. 3 (the concentric lines 
indicate their amplitude in a linear scale and in arbi- 
trary units, say, for instance from at the center to 8 




FIG. 3: Angular radiative cross section (eq.(ll)) and surface 
plasmon cross section (eq.(21) represented by the almost hor- 
izontal lines at a ~ and a ~ 180°), for square defects with 
lOnm side, illuminated by a SPP on silver at 500nm. The 
scale is linear but the units are arbitrary. Each concentric 
line indicates an equal increment of the cross-sections, from 
the minimum at the the center to the maximum at the out- 
ermost. 




FIG. 4: Angular radiative cross section (eq.(ll)) and surface 
plasmon cross section (eq.(21) represented by the almost hor- 
izontal lines at a ~ and a ~ 180°), for rectangular defects 
with lOnm height and 50nm width, illuminated by a SPP on 
silver at 700nm. The scale is linear but the units are arbi- 
trary. Each concentric line indicates an equal increment of 
the cross-sections, from the minimum at the the center to the 
maximum at the outermost. 



at the outermost). For both ridges and grooves er+ and 
<7~ are roughly equal. However in the case of ridges, 
<7p is greater than the maximum value of the scatter- 
ing cross-section into radiative modes da^Jja)/da of 
eq.(ll), by a factor slightly greater than 2. For grooves, 

Op is greater than the maximum value of dcr rad (a) / da 
by a factor slightly smaller than 2. 

Let us now keep the defects height at h = lOnm and en- 
large the width w. Fig. 4 renders the radiation pattern 
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FIG. 5: Angular radiative cross section (eq.(f f )) and surface 
plasmon cross section (eq.(21) represented by the almost hor- 
izontal lines at a ~ and a ~ f80°), for rectangular defects 
with lOnm height and 300nm width, illuminated by a SPP 
on silver at 700nm. The scale is linear but the units are ar- 
bitrary. Each concentric line indicates an equal increment of 
the cross-sections, from the minimum at the the center to the 
maximum at the outermost. 



for a rectangular defect of width 50nm (h = lOnm). The 
emergence of directivity in the out of plane radiation, is 
part of a transitional behavior, in which the radiation 
patterns tend to align and, simultaneously, one of the 
lobes is shrunk while the other is blown up in the ridge 
radiation. Notice that the scattered energy into SPPs 
exhibits the same directivity, going mainly in reflection. 
Eventually, if we keep enlarging the defects until they 
are considerably asymmetric the radiation patterns for 
both ridges and grooves tend to be single overlapping 
lobes (see Fig. 5). Noticeably, the scattering into SPPs 
is greatly reduced. Such similarity is explainable in the 
approximate framework presented in Ref.[22] which turns 
out to be quite acceptable in this limit of large enough 
defects, as we shall substantiate in Sec. VII. In Sec. VIII 
we shall account qualitatively for the transition observed 
in Fig. 4, explaining why the radiation pattern changes 
when the defects are enlarged. 



A. Scattering by square ridges and grooves in the 
Rayleigh limit 

The equivalence between non-elongated subwavelength 
defects and point dipoles gives us a chance to investigate 
in depth the individual radiation pattern of a single scat- 
tering element. 

Fig. 6 shows the averaged the field inside the lOnm ridges 
and grooves, as prescribed in eq.(28) and eq.(29). The 
field induced in a groove is mainly longitudinal while the 
field inside the ridge is mainly transversal. This is due to 
both the illumination and the polarizabilty of the scat- 
terers. When defects are almost symmetric their polar- 
izibilities j3i are nearly isotropic and so the induced field 
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and the incident field are virtually parallel. Hence the 
field induced in a ridge and a groove are nearly paral- 
lel to the incident surface plasmon e spp , which is mainly 
perpendicular to the plane in the vacuum semi-space and 
is mainly parallel to the plane in the metal semi-space. 
Therefore, in the Rayleigh limit, a ridge scatters SPPs 
into radiative modes like a vertical dipole on the plane, 
while the groove scatters them into radiative modes like 
a horizontal dipole on a plane. The results for grooves is 
in agreement with Ref.[20]. 

Interestingly, we also have found numerically in Fig. 6 
that: 



\E 9 J 



E z 



(32) 



especially at short wavelengths. We have devised a vir- 
tual source, that can condense the orientation of the 
equivalent dipole representing a non-elongated symmet- 
ric ridges and grooves. This virtual dipole is defined as: 
q(0) = u x \f\e\ cos 9 + u z sin 6. The fields inside a groove 
and a ridge, are respectively, represented as: 



IE 9 1 



IE' 



|A^q(0)| 



|A £ ^q(90°)| 



(33) 



(34) 



at least as long as eq.(32) holds. 

In reality we can see what happens by means of eq.(24). 
Despite the fact that this equation is only exact for a cir- 
cle in a homogenous background (as explained) we can 
use it to show qualitatively the relation between the field 
inside the groove and the ridge, when their shapes are 
symmetric. If we approximate the polarizability of a 
ridge for that of a circle in vacuum (whose polarizabil- 
ity is calculated through cq.(24)), so 0i = 2/(e + 1). If 
we also approximate the groove polarizability by that 
of a hole in a homogenous metal medium, we have: 
02 = 2e/(e + 1). Hence the field induced inside each 
object is: 



E 



j3ie spp (x = 0,z 

*spp{ 







(35) 



E 5 « f3 2 e spp (x = 0, z = Q-) = f3 2 ^ (36) 



Since these polarizabilities also have the property: 2 = 
efii (the polarizability of a hole in a material is e times 
larger than the polarizability of a particle of the same 
material and the same shape) then \y/e E x \ ~ 
The symmetry of the polarizations and the property 
02 — £0i are strictly true for circular defects in homoge- 
neous media. Our numerical calculations of Fig. 6 shows 
that, even though the field inside a ridge and a groove are 
quantitatively different from those of circular defects in 
homogenous media, the assumption that their mutual re- 
lation is preserved is in very good agreement with the ex- 
act result. Because of the symmetry of the square shape, 
the averaged field inside the square is very nearly parallel 
to the incident field. 
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FIG. 6: The averaged field components (as defined in eq.(28), 
eq.(29)) for a square groove and a square ridge of lOnm side 
in silver, as a function of the wavelength. The scale is loga- 
rithmic with arbitrary units. 



Reflection of surface plasmons square shallow 
defects 



As a corollary of the properties of the fields in a ridge 
and a groove ~ \E X \ we can also substantiate 

that their reflection of surface plasmons is quite similar. 
In fact, we obtain: 



e p 2\^\K 



\E 



(2) | 



p±l 



El 



\E' 



\E 



(i)i 
p±\ 



(37) 



(38) 



Notice that these define a p through cq.(21). Once a ra d 
from eq.(12) and a p are determined the value of the 
transmission of the surface plasmon is a constrained vari- 
able: T = 1 — a~ — <T rac i 7 at least for the lossless case 22 . 
Since a ra d is greater for grooves than for ridges, the 
groove transmission is smaller. 



VI. RADIATION PATTERNS FOR 
HORIZONTAL AND VERTICAL POINT 
DIPOLES ON A REAL METAL INTERFACE 

The first part of the expression eq.(ll) is a pre- factor 
g 3 /S spp whereas the second part is the the radiation pat- 
tern of a point dipole: 



|6W| 2 = *M(a)-p 



(39) 



A groove emits like a horizontal dipole. The angular 
amplitude of the field radiated by a horizontal unit dipole 
p = u x , placed close to the interface z — 0, is & x (a), and 



FIG. 7: Radiative angular intensity \& x \ 2 of a horizontal point 
dipole at an air-metal interface. The radiation patterns varies 
as the metal dielectric constant is varied. The scale is loga- 
rithmic with arbitrary units. 



it does not matter on which side of the interface it is 
placed. $ K (a) can be derived using the relations in the 
Appendix B 2 and the explicit result is: 



2Ve 



cos- 1 a sin a 



ye — cos 2 a + e sin a 



(40) 



and the radiation pattern is |$ K (a)| 2 . Notice $> x (a) 
presents a mirror symmetry about the angle a = 90°, 
the normal to to the plane. Furthermore since <f> x (a) 
never changes sign between and 180° (nor goes to zero), 
the field of a horizontal dipole has one single symmetric 
lobe, where the field always has the same sign. The field 
intensity |$2;(a)| 2 of such lobe is rendered in Fig. 7 for 
different dielectric constants. This radiation pattern of a 
groove shown in Fig. 7, is in agreement with the one rep- 
resented by Ref. [19], obtained with a different numerical 
method. Notice that for lei >> 1: 



2e 



-1/2 



(41) 



That is, when s increases this radiation pattern tends 
to become simultaneously isotropic and vanishing. In 
fact a horizontal dipole does not radiate on a perfect 
conductor 37 . On a small digression it is interesting to no- 
tice an apparent contradiction between treatments such 
as Ref. [38], which considered that a defect in a perfect 
metal were equivalent to a magnetic dipole, while another 
work 20 explains a defect in a real metal corresponds to 
an electric dipole. Actually we have just reconciled the 
two results. We know that a horizontal dipole on a plane 
tends to emit isotropically for large e. This means that 
on a first order expansion in 1/e, the radiation pattern 
of a horizontal dipole on a plane and that of a magnetic 
dipole in vacuum, are identical. 

For finite e the field $ x (a) of a horizontal dipole within a 
real metal would not be thoroughly screened, and while 
the pattern remains symmetric, its isotropy is disrupted 



FIG. 8: Radiation pattern |0 r ad(«)| from the virtual dipole 
q{9) at 8 — 0, 45°, 90°. The scale is linear with arbitrary units. 



parallel to the surface (i.e. a = 0, 180°) to accommodate 
the emergence of the surface plasmons density of states. 
For an individual vertical dipole p = u z , which repre- 
sents a ridge, the angular amplitude of the field is (see 
Appendix B 2) : 



2\e\ s'ma 



ye — cos 2 a + e sin a 



cos a 



(42) 



The field from a vertical dipole also goes to zero at 
a = 0, 180° for a finite e, but since dipoles only radi- 
ate transversally, the field has a third zero at 90°. The 
field is antisymmetric with respect to the normal of the 
plane, while the intensity |$ z (a)| 2 is symmetric, and is 
made up of the two lobes separated by a zero at 90°, see 
Fig. 8. Yet it is important to keep in mind that the field 
of one lobe is in anti-phase with the field of the other. 
Unlike a horizontal dipole, the vertical dipole radiative 
field does not vanish for |e| >> 1 in fact: 



^i 1 ) -> 2 cos a 



(43) 



The total radiation from a vertical dipole has a larger 
weight than the radiation by a horizontal one, by a factor 
of y/e. This can be seen, in fact, from eq.(40) if we assume 
e is large, we get the following relation: 



®£\a) ~ y/e $ x (Q)cosa 



(44) 



In Fig. 8 we represent radiation pattern of q(0) for the 
horizontal and vertical orientations respectively, 9 = 0, 
6 = 90°, which corresponds to our analytic analog of 
the emission pattern of square ridges and grooves respec- 
tively. While we will consider an intermediate orientation 
in the next section, we want to remark here that, due 
to eq.(44), the radiation by both the horizontal moment 
q(0) and a vertical moment q(90°) vanish parallel to the 
plane at a = 0, 180° in a similar manner, as illustrated 
in Fig. 3 

At the same time the far-field emissions of ridges and 
grooves become increasingly different as we approach the 
direction normal to the plane. 
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VII. SOLUTIONS FOR LONG AND SHALLOW 
RIDGES AND GROOVES 

For shallow and long defects w > h and h/X « 1 we 
define the following height-averaged polarization densi- 
ties and fields: 

r h 

P (1 V) = Ae dz'B^(x',z') = 
Jo 

= AeEW(iO h (45) 

where the last equation defines E^^x'). Likewise for a 
groove we can define P( 2 \x') and E^(x') through the 
following equation : 



P (2 \x') = -Ae / dz'E^ix' , z') 

J-h 

= -Ae& 2) (x') h 



-g\*'\ 



(46) 



Notice for |e| >> 1 we can make the approximation 

The benefit of using P^^(x') is that the scattered-field 
coefficients for these defects in the far zone, B^ d (a) and 
@p± ' are those emitted by a chain of point-dipoles on 
the surface over the segment w, and set at + and _ for 
ridges and grooves, respectively. 



The scattered field angular amplitude Q^ad ( a ) 
eq.(10) and eq.(13) is obtained as: 



from 



Jo 



■xgx cos a 



(47) 



This holds for the scattering into surface plasmon modes 
as well since we have: 



u p± — 



(0) 



dx' e Tik ** x ' P {v \x'). (48) 



When we illuminate a shallow and long defect, with a 
SPP, an equivalent linear density of dipole sources P(x') 
stems from how the induced fields are distorted inside the 
scatterer, namely by its polarizability. When the defect 
is larger in the horizontal direction than in the verti- 
cal one, ridges and grooves were found to give the same 
scattering by an approximated Rayleigh expansion 22 . We 
have an alternative first principles argument to justify the 
Rayleigh expansion result, which is based entirely on the 
assumption that these defects are needle shaped. The 
field induced in these defects tends to be that induced 
in a needle-shaped protrusion placed horizontally on the 
surface + in the case of a ridge. For a groove we have a 
horizontal needle-shaped cavity at 0~. In such idealistic 
simplification it is clear-cut to deduce the fields inside the 
defects from the boundary conditions. Namely the par- 
allel component of the incident field is always continuous 
and equal, as in eq.(Al)and eq.(A2): 
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|4 CI) W| 
|e|K"'(*)| 

Groove 

" " " |£ Z °'M| 




100 150 

x(nm) 
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E lx (x')=e sp (x',0) 



E 2x {x') 



(49) 



FIG. 9: The averaged fields component inside of a ridge and 
a groove, \ the quantity eE^ and the rest of the compo- 
nents Ex , Ei , E^\ for rectangular defects of w = 300nm 
and h = Wnm. The system is illuminated by a SPP in loss- 
less silver at A = 700nm. The scale is linear with arbitrary 
units. 



which preserves the continuity of eq.(14). However, we 
are generating fields which, normal to the surface, make 
up for the discontinuity perpendicular to the metal sur- 
face of eq.(15). In fact, for a horizontal needle-like ridge, 
the boundary conditions imposed by the continuity of the 
displacement vector are: 

E u {x') = e spp {x', 0+) • u z /e = l/e (50) 

while for a needle-like slit: 

E 2z (x') = ee 3pp (x' > 0-)-u z = l. (51) 

Ultimately: 

E lx {x') = E 2x (x') (52) 

sE lz {x') = E 2z {x') (53) 

which, matched with eq.(14) and eq.(15), yields: 

G(R,x',z' = 0+)-Ei(x') = (54) 

G(R,x',z' = (T)-E 2 (x') 

and thus the property of producing the same scatter- 
ing coefficients, previously found in Ref.[22]. Of course 
this is just an approximation, but it explains why elon- 
gated defects have similar scattering properties. In real 
life the plasmon scattering by protrusions and indenta- 
tions is similar because, far from the edges, a shallow but 
elongated defect behaves as an infinitely elongated one, 
as confirmed by numerical calculations. As an example 
we report in Fig. 9 a numerical calculation of the fields 
averaged over the height for defects of w — 300nm and 
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h = lOnm. This shows that eq.(52) and eq.(53) are quite 
accurate at the center of the defect, and deviate from 
the needle model prediction due to fringe effects at the 
edges. 

It is worth mentioning that this equivalence is valid in 
the Rayleigh limit when the defect size is much smaller 
than the wavelength, and may be altered at resonant 
wavelengths. 



VIII. THE TRANSITION FROM SHORT AND 
SHALLOW DEFECTS TO LONG AND SHALLOW 
DEFECTS: OBLIQUE DIPOLES ON A REAL 
METAL PLANE 

Everything we just said for symmetric surface defects 
was based on the fact that their aspect ratio equals one. 
As the defect width is increased, the aspect ratio be- 
comes larger and this leads, progressively, to an asym- 
metric polarizability tensor. The first effect is that the 
field induced is gradually less and less parallel to the 
incident field. Therefore a ridge would develop a non- 
negligible horizontal electric field component, thus ceas- 
ing to be equivalent to a vertical dipole. Likewise the 
groove, which in the symmetric case behaves as a hor- 
izontal dipole, gradually starts having a non-negligible 
vertical component as its shape is elongated. The pro- 
cess goes on until we recover the case of a needle shaped 
defect of section VII. The fields inside a defect having in- 
termediate width, as in Fig. 4, are intermediate between 
those for the needle case and the square symmetric case. 
Therefore in these cases defects emit qualitatively like 
oblique dipoles, with the orthogonal components out of 
phase. 

In order to understand better the radiation pattern by 
ridges and grooves we decompose the oblique dipole in 
its horizontal and vertical components. 
First of all, we focus on the mechanisms involved radi- 
ation pattern for a ridge v = 1. From eq.(30) a dipole 
with arbitrary orientation emits close to the surface, with 
a field angular amplitude: 

©™» = * (1) («)-P (1) = <M«)A«(a) (55) 

where A^ (a) — + (^i 1 \a)/^ x (a)j p^ and equals: 

A^(a)= P ^ +l £C 2 ° Sa P U (56) 
ycos^ a + \e\ 

A^(a) shows that the contribution to the radiative field 
coming from the vertical and horizontal dipole on a metal 
plane have a phase difference of 90°. This was already 
evident from eq.(44), when e < 0. Such phase dif- 
ference arises from the impedance of a metal plane 22 

z s = -i/yT4 

The radiation pattern for a dipole with arbitrary orien- 
tation and lying above the metal, is written in our for- 
malism as: |$ a; (a)A( 1 )(a)| 2 . The net angular amplitude 



p x =expM) 




FIG. 10: Radiation pattern |6 r ad(«)| for a point dipole: p = 
u^e"^ + u z , lying on top of a metal surface. The scale is linear 
with arbitrary units. 



for an oblique dipole is resolved into the superposition of 
the angular envelope of the horizontal dipole (shown in 
Fig. 7), with the other radiation factor |AW(a)| 2 . This 
last factor contains both the orientation and phase of the 
field. To envisage how these combine we may develop 
lA^^a)! 2 into three terms. These consist in the individ- 
ual emission from the horizontal and vertical dipole plus 
an interference term: 

|A (1) (a)| 2 = b.| 2 + ^^b,| 2 + 
|e| + cos z a 

-2 | £| 3 [ PlxP * z ] cos a (57) 
^/cos^ a + \e\ 

In the presence of the plane metal background, we have 
that horizontal and vertical dipoles behave as individual 
sources but their interaction presents an intrinsic added 
phase difference of 90°, which is due to the different in- 
teraction of a horizontal and a vertical dipole with the 
plane. As a result, when in phase they do not interfere, 
and their radiation pattern is always symmetric regard- 
less of the orientation of the dipole. This is the case for 
q(45°) where, as in Fig. 8, the radiation pattern is the 
sum of the angular intensity of a vertical and a horizontal 
dipole, so that at 90° there is a minimum due to the van- 
ishing of the vertical dipole contribution, and yet never 
goes to zero because of the horizontal dipole contribu- 
tion. Nevertheless, when the dipole components are not 
in phase, we can get asymmetric radiation patterns and 
additional zeros (to those at 0° and 180°), because the 
interaction term can be negative. In such case the in- 
teraction of the horizontal radiative field (with only one 
lobe) with the vertical radiative (with two lobes of dif- 
ferent sign) is responsible for an asymmetric radiation 
pattern and exhibits directionality. This is illustrated in 
Fig. 10 for a dipole emission whose main contribution 
comes from the vertical dipole. In Fig. 11 we show the 
radiation pattern for a dipole whose main contribution 
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p. =8 




FIG. 11: Radiation pattern |O ra( j(cf)| for a point dipole: p = 
eu x + e ll ^u z , lying on top of a metal surface. The scale is 
linear with arbitrary units. 



comes from the horizontal dipole radiation. 

For the case of a grooves [y = 2), the radiative angular 

field amplitude is, from eq.(17): 



IX. CONCLUSIONS 

Our analysis of the surface plasmon scattering by 
square shallow defects into radiative modes and plasmon 
modes, reveals that a groove scatters more of the incident 
energy than a ridge does. The reflection by a symmetric 
ridge and a groove is similar and so is the radiative emis- 
sion close to the horizontal direction. Indeed their scat- 
tering essentially differs in the vertical direction, where 
a groove scatterers while a ridge does not. When defects 
start to become longer in width we saw the polarizabil- 
ity gets more asymmetric. Correspondingly, since both 
components of the incident plasmon are out of phase, de- 
fects are equivalent to interfering horizontal and vertical 
dipoles on a plane, which interfere constructively in some 
direction, thus producing directionality in the radiation 
pattern. Finally when ridges and grooves are shallow and 
long they tend to produce the same scattering as, apart 
for fringe effects, their polarizability exactly counterbal- 
ances the discontinuity of the incident surface plasmon 
field at the air-metal interface. 
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ei.a» = <MaW 2 )(a) e -^'l^ (58) 
\/cos z a + \e\ 



where remember we have also added the approximation: 
kf ~ gyfi for |e| » 1. 

Remarkably, as opposed to the the dipole emission over 
the surface, in the net emission from a dipole under the 
surface the horizontal dipole contribution has a greater 
weight than the vertical dipole contribution. Apart from 
this, all the arguments used for a dipole over the surface 
apply. 

The interaction between the vertical and horizontal com- 
ponents of the field induced in the field generates the 
directional patterns of Fig. 4. For a ridge with length 
slightly larger than its height the directional radiation is 
dominated by its vertical component. Fig. 10 exempli- 
fies the effect of the interference of a dominant vertical 
component with a smaller but non-negligible horizontal 
component. For even larger aspect ratios the contribu- 
tion from the other component may be comparable. 
Likewise when a groove has a small aspect ratio it is pre- 
dominantly a horizontal source interfering with a smaller 
vertical source. The result is in an interference pattern 
that looks like the one rendered in Fig. 11. Yet again 
this can be modified by increasing the aspect ratio. This 
transition is in good agreement with Fig. 11 of Ref.[19] 
where, using a different numerical method, the radiation 
pattern of a groove was computed for different aspect 
ratios. 
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Appendix A: Surface Plasmon Polariton Mode 

The incident illumination is the field of a surface plas- 
mon wave mode propagating in the positive x direction 
(+) or negative x direction (— ) is: 

e%;l\r)=(^- + u^j e^M+V*), z > 0j (A1) 

e£2>(r) = + e *(±^+CW), z < 0. (A2) 

where k p = g(s/(e + l)) 1 ^ 2 , k pz = ig/y/—s — 1 and 
k™ = —ekp Z . This can, alternatively, be written as: 

e spp ±{r) = for z > 0; and e spp± (r) = e ( spp %) for 

z < 0. 

The magnetic field associated is continuous at the inter- 
face and equal to: 

— i 

h spp ±(r) = — V x e spp± (r) (A3) 

Now consider a lossless metal, characterized by a real 
and negative dielectric constant e and consider a plasmon 
moving in the forward direction, (the subscript + will be 
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omitted). The incident Poyinting vector of the plasmon 
in the air side is: 



w=l _ / 
spp / 

Jo 



az e spp x ti spp ■ u x = — — — - 

9 z l' t ; 



1 



(A4) 



pz 



while in the metal is S" pp = f ° dz e spp x h* • u x — 

ZfSg^p 1 where Z s = —ij \f\e\. The total Poynting vector 
energy flux associated to a plasmon mode in a lossless 
metal is: 

c _ ou=i , q u=2 _ (e + l)(e 2 - 1) 
Jspp — >J spp t o spp ~ 2g £ 3 



Appendix B: P-Modes 



2. Key Identities 

The following expressions for the reflection and trans- 
mission coefficients are essential to derive cq.(40) and 
eq.(42): 



4 1,2 V) 

1 + r p (a) 
1 - r p (a) 





2y / esina 




- cos 2 a + e sin a 




2Ve — cos 2 a 


Ve 


— cos 2 a + e sin a 




2e sin a 




— cos 2 a + e sin a 



(B5) 

(B6) 
(B7) 



We shall repeat, out of completeness, the explicit ex- 
pression for p- waves, particularly in the far field when 
k/g = u r . In this case these modes are expressed in 
terms of the direct space polar angle a by noticing that 
k x = g cos a and k z = ki"~^ = gsina in the air semi- 
space and fc™ = k z u ~ 2 " > = g\Je — cos a 2 in the metal. 
Hence 

kp(a) = i (k z u x =F k x u z ) = sinau^ =p cosatLj 



(a) = 



{k™u x T k x u z ) 




cos a 
u x T — — u z 



(Bl) 



1. Reflection and Transmission coefficients for a 
plane surface 

For reference, we give here the Fresnel coefficients for 
an air metal interface. In the present treatment we only 
deal with the reflection coefficient for a p-wave propagat- 
ing from air to metal, and this is : 



r _ r (i,i) 



ek z 



ek z 



(B2) 



where notice that, for the sake of tidiness, we omit the 
superscript throughout. 

As to the transmission coefficients the one for a wave 

(2 1) 

(2,1) propagating from the metal to air is t p ' , while 
the one for a p-wave transmitted from the air medium to 
the metal is tp 1 ' 2 ^. 

,(2,1) _ 2k™T/e ( 1)2 ) _ 2k z ^/e 

P Urn _i_ p U V Urn i F U • ' 

rv z T c r\j Z fv z c rv z 

Notice that the transmission coefficients are related as 
follows: 



,(1.2) ,(2,1) 
bp bp 



.(1) 



<„( 2 ) 



(B4) 



Appendix C: Asymptotic Green's Tensors 

The asymptotic expressions for the Green tensor for 
3D scatterers are found in references 11,33,39 . We have al- 
ready presented the derivation scheme for bi-dimensional 
defects in Appendix B of Ref. [30] , for a groove. As ex- 
plained therein the Surface plasmon Green tensor and 
the far-field Green tensor are obtained from its angu- 
lar spectrum. From the relevant Sommerfcld integral the 
surface plasmon contribution is obtained by applying the 
residue theorem and the far-field Green tensor instead is 
obtained by applying the method of the steepest descent. 
For the case of the ridge we use the total Green tensor of 
the background in the vacuum semi-space. This can be 
written as the sum of the direct Green Tensor (the free 
space green tensor ) and the indirect green tensor (which 
gives the contribution due to the reflections at the metal 
plane interface). Hence 



G« (R, r') = G (R, r') + G S (R, r') 



(CI) 



where the spectral representation for the direct Green 
tensor is: 

Go (R,r') = A |°° ^ e ^(W W W) k + k + 

(C2) 

while for the indirect Green Tensor: 

G s (R,r') = ^ J°° ^e^W) eiM z + .') rpk - k + 

(C3) 

Applying the residue theorem and the steepest descent 
method to G^'(R, r') we end up with eq.(4) and eq.(18) 
for (y = 1). 

For the groove case we need to expand the Green Tensor 
connecting a point in the metal to a point in air. This is 
just: 

G ( 2 )(R,r') = — [°° % e ik x (x- x ') e i(k z z-k?z>) x 



,(2,1) k ™+ k + 

L p K p *-p 



(C4) 
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Applying the residue theorem and the steepest descent 
method to G^(R, r') we end up with eq.(4) and eq.(18) 
for (u = 2) . Notice that the form of G°°(a,r') = 
— <I>Wu Q given in Sec. Ill, is obtained by recognizing 
kp = u Q . 

One more subtlety, that might be confusing, is how we 

(2 1) 

pass from the transmission coefficient t p in the inte- 

(1 2) 

gral to the transmission coefficient tp ' in the asymp- 
totic form 3>( 2 ) . This comes about because when we 



apply the method of the steepest descent to the integral 
we get eq.(4 ) with: 



qoc ( r ') = ^ k ™ Ua = t (i,2) k - Uq (C5) 



where, in the last equation, we have used the identity 
eq.(B4). 
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